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Abstract 

We studied the current-induced magnetization dynamics of a domain wall confined in a magnetic 
wire with bi-axial anisotropy. We showed that above the threshold current density, breathing-mode 
excitation, where the thickness of the domain wall oscillates, is induced by spin-transfer torque. 
We found that the breathing-mode can be applied as a source of microwave oscillation because 
the resistance of the domain wall is a function of the domain wall thickness. In a current sweep 
simulation, the frequency of the breathing-mode exhibits hysteresis because of the confinement. 
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Recent advances in spin electronics have revealed that the current flowing through a 
magnetic nanostructure with a non-collinear magnetization configuration can excite magne- 
tization dynamics. Current-induced microwave generation has attracted a lot of attention 
because it will be a candidate for applications in future wireless telecommunication tech- 
nologies. Most studies of current- induced microwave generation have been carried out in 
magnetic multilayersi^ 1 ^^^. In these experiments, uniform precession of the free layer 
magnetization*^ is driven by spin-transfer torque^^*, and the motion of the free layer 
magnetization is measured using the CPP-GMR or TMR effect. 

Only a few works have theoretically suggested on the current-induced microwave genera- 
tion of a domain wall^ 1 ^. He and Zhang proposed an application of the oscillating motion 
of the domain wall under current and magnetic field as a source of microwave oscillation^. 
Ono and Nakatani proposed a microwave oscillator using the rotating motion of the domain 
wall^. 

On the other hand, it is known that a domain wall produces an additional resistance in 
magnetic wires. According to the theory of Levy and Zhang*^, the resistivity of a domain 
wall is inversely proportional to the square of the domain wall thickness. Therefore, if we 
excite a breathing-mode where the thickness of the domain wall oscillates by application of a 
dc current, we can use the oscillation as a microwave source. From the view point of physics, 
it is also important to study the current-induced magnetization dynamics (CIMD) of the 
geometrically confined domain wall since CIMD is in general different from magnetic-field- 
induced magnetization dynamics and little is known about the CIMD of a geometrically 
confined domain wall. 

In this paper, we investigated the CIMD of a domain wall confined in a magnetic wire 
with bi-axial anisotropy by a confining potential due to the wire shape shown in Fig. QJa). 
We showed for a confined domain wall that above the threshold current density, spin-transfer 
torque induces breathing-mode excitation*^ 4 , where the thickness of the domain wall os- 
cillates. The current-induced breathing-mode can be applied for a microwave oscillation 
because the resistance of the domain wall is a function of the domain wall thickness. We also 
found that if the current density is adiabatically changed, the frequency of the breathing- 
mode shows a hysteresis loop because the confining potential enables the breathing and 
pinning states to coexist below the threshold current density. 

The system we consider is a 180° domain wall confined by a potential due to the wire 
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XIX 

FIG. 1: (a) Bloch and Neel walls confined in a magnetic wire are schematically shown. The 
regions outside the wire are represented by shading. The solid and hollow arrows represent the 
magnetization vector of the Bloch and Neel walls, respectively. The dotted circles represent the 
uniform rotation where the exchange energy of the domain wall takes a constant value, (b) Potential 
energies are plotted against the width of the domain wall A for (fro = (solid line) and (fro = tt/2 
(dashed line). The horizontal axis is normalized by the width of the domain wall of the ground 
state for (fro = 0. The energy minima are indicated by arrows. 

shape, as shown in Fig. [T](a). For simplicity, we modeled the system as a one-dimensional 
domain wall along the x-axis with a confining potential, where the directions of the magne- 
tization vectors are represented as M = (sin 9 cos 0, sin 9 sin 0, cos 9), where 9 and denote, 
respectively, polar and azimuthal angles for each position. We assume that the system has 
a bi-axial anisotropy such that the z- and y-axes are easy and hard axes, respectively. 

Let us begin with a brief introduction of the theory of a one-dimensional domain 
waiy^^L^. i n the absence of current and confining potential, the energy of the system is 
expressed as 

/oo 
dxJ [9' 2 + sin 2 #0' 2 ] + fT c sm 2 9 [l+r sin 2 0] , (1) 
-oo 

where the primes denote differentiation with respect to spatial x coordinate. The first term 
is the exchange stiffness energy with a stiffness constant of J. The second term represents 
the bi-axial anisotropy characterized by the anisotropy constant for the easy axis, K e , and 
the ratio, r = Kh/K c , where is the anisotropy constant for the hard axis. 
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It should be noted that if the system has uni-axial anisotropy, i. e., r = 0, the ground 
state is degenerate with respect to the azimuthal angle (ft which is independent of the spatial 
coordinate. Since the magnetization vectors of the top and bottom electrodes are aligned 
parallel to the z-axis, the Bloch and Neel walls correspond to the azimuthal angle of (ft = it/ 2 
and = 0, respectively. The magnetization vectors of the other ground states lie on the 
dotted circles as shown in Fig. []Ja). The polar angle configuration 9 of the ground state is 
given by 



x 



6(x) = arccos tanh J , (2) 

where A = J/K c is the domain wall thickness. The degeneracy of the ground state with 
respect to the azimuthal angle (ft is broken if the system has bi-axial anisotropy, i.e., r ^ 0. 
Assuming that the azimuthal angle (ft is independent of the spatial coordinate, the polar 
angle configuration of the minimum energy state for each value of (ft is expressed by Eq. (j5J) 
with A = A (l + rsin 2 0) X ^ 2 , where A = a/ Jj K e represents the domain wall thickness of 
the ground state with = 0. 

Since we are not interested in the domain wall propagation along the wire but in the 
breathing-mode excitation, we assume that the domain wall is confined in a certain region 
of the wire by a confining potential V c {. We also assume that the characteristic length of 
V c f is much larger than A. We adopt <fto, X and A as collective coordinates^ 2 ^. Here (fto is 
defined by 

dx(ftsm 2 9, (3) 



and X is defined through the polar angle configuration of the ground state, 9o, given by 

'x-X 



9n = arccos 



tanh 



(4) 



A 

where X denotes the position of the domain wall center. Substituting Eqs. (3) and (4) into 
Eq. (1) as 9 = 9q + (9 — 9q) and <ft = <fto + (<ft — <fto), and ignoring the spin-wave excitation of 
9 — 9 and (ft — (fto, the energy of the domain wall is obtained as 



2 7 

Efa, A) = -j- 



j + j (l +rsin 2 o ) 



(5) 



In Fig. d](b) we plot the domain wall energy of Eq. (jSJ) with r = 0.3 as a function of the 
normalized domain wall thickness A/ A. The solid and dotted lines correspond to (fto = and 
7r/2, respectively. As shown in Fig. QJb), the value of A which minimizes the domain wall 
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energy depends on 0o- Therefore, if the precession of the domain wall around the azimuthal 
axis is induced by an applied current, the oscillation of A and therefore the resistance of the 
domain wall are also induced. 

For simplicity we assume that the confining potential takes the form V c f = 
F L [sinh 2 (X/L)] , where F and L represent the magnitude and the characteristic length of 
the confining potential. As mentioned above we also assume that the characteristic length L 
of the confining potential is much longer than the thickness of the domain wall A. Thus the 
domain wall thickness is not related to L and the domain wall is different from the so-called 
geometrically confined domain wall whose thickness is determined by Lr^. 

In order to systematically derive the equation of motion described by the collective co- 
ordinates, we adopt the Lagrangian met ho cl2ii22i2&i21 The Lagrangian corresponding to the 
torque exerted on magnetizations of the domain wall by an applied current^ 1 ^ is given by 

-/>(^^i^' [1 - cos, ' , ■ < 6 > 

were £ is the ratio between the precession time due to the exchange interaction and the spin 
relaxation time for spin accumulation; e is the electric charge of an electron; j e , the charge 
current density; P, the spin polarization of the charge current and 7, the gyromagnetic 
constant. Then the total Lagrangian of the system under the applied current 2 ^ 2 -* 2 ^ is given 
by 

1 f 00 

£ — — dx(4> — j(f)') [1 — cos 6] 
1 J —00 

/oo 
dx [9' 2 + sin 2 9 <f)' 2 ] 
-00 

/oo 
dx sin 2 9 [l + r sin 2 <p] - V ci , (7) 
-00 

where j = /iB-Pj c /(l + £ 2 )e represents the spin-current density. The dots denote differenti- 
ation with respect to time t. In order to obtain the effective Lagrangian described by the 
collective coordinates X, <fro and A, we substitute 9q and <fro for 9 and in Eq. (J7J), respec- 
tively, and then perform integration with respect to x. We obtain the following effective 
Lagrangian, 

2 • 2 7 

C = -(X+ j)c/> - — - 2XK d (1 + r sin 2 O ) - V*. (8) 

7 A 

In the Lagrangian formalism, the effect of the Gilbert damping 2 ^ is conventionally taken 
into account by the Rayleigh dissipative function method^ 2 - 1 ^. The dissipation function is 
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defined by 



a f j \/n P -n>\2 i • 2/1/1 -J.'\2 



T=-Jdx^{9- <- 2 ey + sin 2 9{<j> - } , 



where a is the Gilbert damping constant. The terms proportional to (3 reproduces so- 
called /J-term^k&is&i&M and describes drift effect^. One can easily confirm that Eq. ([9]) 
reproduces the torques coming from the Gilbert damping and (3 terms as —M x SJ-'/SM = 
—aM x M + (3jM x M' . In order to obtain the effective dissipation function described 
by the collective coordinates X, <p and A, in the same way as the derivation of the effective 
Lagrangian, we substitute 9q and <po for 9 and 0, respectively, and then perform integration 
for the coordinate x. We obtain the following effective dissipation function, 

r a \X 2 - 2 vr 2 A 2 2/5 • 1 , , 

J r=_J T+ ^ + __. + _£,*j, (10) 

where the terms which is independent of the time derivatives of the collective coordinates 
are dropped. 

The equation of motion for the domain wall is obtained by using the effective Lagrangian 
of Eq. (jSJ), the effective dissipation function of Eq. ffTOl) and the Euler-Lagrange equation, 

d5C_5C 5F 

di~5q~ ~ Jq~ ~~ &p ^ ' 

where q being in X, <p and A. After some algebra we obtain 

0o + «f = -^j + jF*(X), (12) 



- X + a\<j) j = j - 7^-Asin20 o , (13) 
ii 2 a ■ J \X A / r , , s\l , ,s 

V = i[a-a( 1 + 2( 1 - C082 «)J- < 14 > 

where F c f = — 9xKf- Equations f|T2l) . ffT3l) and (TH1) reduces to Slonczewski's equation 
of domain wall motion^'^I when the dissipation for dynamics of A is neglected. The (3- 
proportional term in Eq. ()12p represents the torque from the so-called /?-term of the Landau- 
Lifshitz equation^. We note that A is a dynamical variable and is independent of 0o in 

Eqs. flED-QH^. 

We performed numerical simulations based on Eqs. ( fT2l - (fT4l) . The equations were solved 
using the implicit Runge-Kutta method. Since the confinement force, F* , is much larger 
than the /3-term in the present situation, we ignored the first term of the right-hand-side 
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FIG. 2: (a) The position of the domain wall X is plotted as a function of time t. The current 
densities are taken to be j =0.37, 0.31, 0.25 and 0.19 for top to bottom. The unit of the current 
density is taken to be 7J/A. (b) The angle is plotted as a function of time t. The plot for 
j =0.19 lies on the horizontal axis. In both panels, the unit of time is taken to be A 2 /(7<7). 

of Eq. fl!2p . We took the parameters of the confining potential to be L = 2. OA, and Fq = 
100.0J/A 2 to efficiently confine the domain wall to the confinement region. This condition 
corresponds to the case in which the cross section of the wire exponentially increases tenfold 
in the linear dimension by the displacement to X = L from 0. We also set the Gilbert 
damping constant at 0.01 to reproduce typical experimental systems. At initial time t = 0, 
the thickness, position, and azimuthal angle were taken to be A = X,X = 0, and 0o = 0, 
respectively. 

In Figs. 12(a) and [2(b) we plot the position X and the angle 0o of the domain wall as a 
function of time for various values of j, respectively. The current is switched on at t = and 
then is kept fixed. As long as the spin-current density j is smaller than the critical value 
j c , the position of domain wall X showed little deviation from its initial value of X = 0. 
In our simulation j c = 0.22 7J/A. Hereafter, the unit of the spin-current density is taken 
to be 7J/A. As shown in Fig. [2(b), the angle 0o also shows little deviation from its initial 
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FIG. 3: Thickness oscillation, (a) The normalized width of the domain wall A/A for j=0.25 is plotted 
as a function of time t. (b) The power spectrum density P{oo) of magnetoresistance obtained by 
using the theory of Levy and Zhangi^. Solid, dotted and dot-dashed lines correspond to the current 
density of 0.25, 0.31, and 0.37, respectively. The unit of the current density is taken to be 7J/A. 

value 0o = for j = 0.19 < j c . Above the critical current, j > j c , X moves to a certain 
position which is determined by the competition between the confining potential and the 
spin-transfer torque. The angle 0o linearly increases with increasing time and the domain 
wall precesses around the x-axis. 

The depinning of the angle 0o shown in Fig. [2](b) is similar to Walker's breakdown^i^i. 
However, the j c we obtained is not equal to Walker's threshold jw = rJ/X. As we shall 
show later, the difference in j c and jvv reflects the essential difference in dynamics between 
confined and unconfined domain walls, such that coexistence of the oscillation and pinning 
states for j < j w depends on the confining potential. 

Since the system has bi-axial anisotropy, the precession of the domain wall induces the 
oscillation of the domain wall thickness A, the breathing mode^, as shown in Fig. [3](a). 
According to Levy and Zhang's theory^, the resistance AR of a domain wall depends on 
its thickness as AR ~ 1/A. Thus, the resistance AR oscillates due to the breathing mode. 
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In Fig. Mjo) we plot the power spectrum density of (A/A) ~ AR defined as 

P M = J d re-^J itm £— y (15) 

For j > j c , the oscillation in AR was observed as expected above. For the typical exper- 
imental situation the frequency is on the order of several tens GHz^. For each value of 
j, P{oj) has a single sharp peak. This means that the system is a useful candidate for a 
microwave source. The peak frequency of the power spectrum is proportional to j, which 
means that we can control the microwave frequency by the current. 

The intensity at the peaks is on the order of 0.01 and depends on the amplitude of the 
thickness oscillation. In this case the amplitude of the thickness oscillation is on the order of 
O.lA. The value of the intensity is consistent with the amplitude of the thickness oscillation 
because P(ui) is proportional to the square of the amplitude normalized by A. The amplitude 
of the oscillation is about Ar. Thus, the intensity is controlled by the ratio of anisotropy 
constants, r, and does not depend on j. We note that A is always shorter than A as shown 
in Fig. [3]^a). Namely A does not agree with the thickness minimizing Eq.Q because of the 
effect of the Gilbert damping in Eq. (j!4p &. 

Next we move onto adiabatic current sweeping. The current was adiabatically increased 
from to 0.4 which is above Walker's threshold current jw, an d then was adiabatically 
decreased from 0.4 to 0. From Eq. (fl~4l) . the frequency of the breathing mode in the steady 
state is related to the velocity of <fto as /o = ^o/ 71 "- Figure H] shows the breathing-mode 
frequency fo- As the current increased from 0, the frequency was kept at zero below the 
threshold current j' c = jw^min/A, where A min = min A = Xy/1 — r. When the current reached 
j = j' c in the simulation, the frequency jumped to a certain value and then linearly increased 
as the current further increased up to 0.4. As the current decreased from 0.4, the frequency 
linearly decreases down to a much lower current than j' c . That is to say, below the threshold 
current, the breathing and pinning states coexist. The behavior differs surprisingly from 
that <^0o^ ~ V ' 3 2 ~~ j\\T^ above the threshold current j w in Walker's theory, where a unique 
state is permitted for each current. 

The reason for current-dependence behavior of the breathing motion frequency is that 
the confining potential enables the breathing and pinning states to coexist. In fact we can 
see two time-averaged solutions of Eqs. (|T2|) - (jl~4l) : the pinning solution is < X >=< 
X >=< <fi >= 0, < 0o >~ arcsin(jA/7r J)/2 and the breathing solution is < X >= 0, 
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FIG. 4: The breathing-mode frequency /o is plotted against the current density j. A hysteresis 
loop appears with increasing and decreasing applied current density (indicated by arrows). 

< 0o >~< j/olX >~ 7 < Fg(X) >. As mentioned before, the threshold j c originates from 
the coexistence of the two states. We note that if the confining potential does not exist the 
breathing mode vanishes because < F^(X) >= 0. The existence of the confining potential 
induces a drastic effect in the domain wall motion. 

In conclusion, we examined the current-induced magnetization dynamics of a domain 
wall confined in a magnetic wire with bi-axial anisotropy. We showed that breathing-mode 
excitation, which produces resistance oscillation, is induced by spin-transfer torque. The 
result means that the confined domain wall is a powerful candidate for a microwave oscillator. 
We also found that the dependence of the frequency of the breathing mode on the current 
shows a characteristic hysteresis loop originating from the confining potential. 
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